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Abstract — The influence of neutron irradiation upon 
(nonlinear physical) elastoplastic deformation of components of 
constructions under cyclic force disturbance is considered. As an 
example the problem of cyclic irradiation-force winding 
consider unsymmetrical with respect to thickness sandwich 
beam with external bearing layers made of metal and 
incompressible with respect to thickness internal layer (filler) 
made of polymer. Well effects resulting in appearance of 
additional volume deformation. Furthermore the differences 
used during the numerical solving reflects of the nature of the 
physical structure of body under influence the flux. On basis of 
experimental data the formula of irradiation reinforcement is 
suggested. The result may be extended in case of any given 
n-cyclic loading The theorem about variable loading is proved. 
Theorem about cyclic loadings of elastoplastic body in neutron 
flux allows to simplify essentially a whole class of boundary 
problems. Numerical solving for sandwich nonlinearly beam is 
adduced. 

Index Terms — irradiation-force, neutron flux, numerical 
solving, , variable loading. 

I. Introduction 

Radiation treatment of rigid body is conducted by 
numerous effects resulting in appearance of additional 
volume deformation, changing of elastic and especially 
plastic properties of substance. Therefore it is necessary to 
bring in appropriate corrections into posing and solving of 
boundary-value problems concerning single and cyclic 
loading of elastoplastic components of constructions. The 
main factors are irradiation reinforcement of substance 
(increasing of yield point) and irradiation swelling (increasing 
of volume deformation). Below an attempt to extend the 
theory of variable loading by Moskvitin [1] upon discussed 
class of boundary-value problems. 

II. Procedure of methods of solution. 

A. Volume Deformation and Irradiation Reinforcement 

Let’s consider initially homogeneous isotropic body, 
occupying half-space Z ^ 0. If upon the border ( Z = 0 ) 
parallel to the axe z are felling neutrons with identical 
average energy and intention (p 0 , the intention of neutrons, 

reaching the plain Z = const , will be [2], [3] (p(z) = (p 0 e - ^ . 
The value ft is called macroscopic effective section and has 
valuation of 1/m. 


If <p Q is not dependent upon time, up to the instance t 
through the section Z will pass the flow 

I{z) = <\>fte vz . (1) 

In rough estimation we may consider, that changing of the 
volume of substance is directly proportional to the flux I(z ) , 
and consequently 0 7 = B/(z) , where B— is experimental 
constant. The value 7 0 = (p {} t will give total flux of neutrons 
upon unit of area of the surface of the body. In reactors (p {) 
has the value of 10 17 -10 ls neutrons/(m 2 s), at that I Q can 

reach the value of 10“' -10“ neutron /m , and u I approaching 
to value 0.1. Consequently, depending upon of energy of 
neutrons and properties of irradiated substance B can has the 
value of 10" 28 — 10" 24 m 2 /neutron. 

Dependence of the modulus of elasticity, yield and solidity 
points and whole tension diagram upon 7 0 for different 

energies was experimentally investigated after irradiating of 
samples in nuclear reactors. Results of experiments show that 
as usual the modulus of elasticity changes weakly (increase by 
1.5-5 % relatively not irradiated sample). As for yield and 
solidity points - they are very sensitive to the irradiating and 
the yield point especially. 

For massive bodies with flat boundary the number of 
neutrons passing at the depth z under this boundary is 
estimated by the formula (1), that's why the yield point will 
vary along z ■ At the surface of the body ( Z = 0 ) the 

influence of irradiation upon the plastic limit <7 y is 

satisfactory characterised by the formula of irradiation 
reinforcement [2]: 

G v =a ro[ 1 + A ( 1_ex P(-^o)) 1/2 J (2 > 

Here cr y0 - the plastic limit of not irradiated substance. At 
the depth Z the formula takes the form 
o y = o y0 [l + A(l - exp (- fl )) 1/2 ] , 

where the value of neutron flow I{z) is characterised by 
formula (1). Let's denote appropriate values of deformation as 
£ y0 ’ £ s ’ - are substance constants, taken from 

experiment. For example, if for aluminium alloy we accept 

A = 1.09 , £ = 9.73x10 ~ 6 m 2 /neutron, then the fig.l 

indicates the satisfaction with known experimental data. Dark 
points - experimental data, solid line - the estimation 
by formula (2). 
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0 2 4 6 /TO 24 , neutron/m 2 

Fig. 1 experimental data, where dark points - experimental 
data, solid line - the estimation by formula (2) 


B. Statement of problems of the theory of small-scale 
elastoplastic deformations 

Let’s consider the process of complex influence of external 
force and radiation loading upon deformable body within the 
theory of small-scale elastoplastic deformations. Suppose, 
that the body in natural state is influenced simultaneously by 
instant external forces F' , R' with boundary moving u' jo 

and neutron flow with the value / 0 = (pt . It is supposed, that 

under such influence areas of elastic and plastic deformations 
appear in the body. Let’s neglect the changing of modulus of 
elasticity caused by the neutron irradiation. Arising within the 
body stresses deformations and moving will be marked by one 
upper dash. 

In elastic areas of solid body Hooke's law is valid and 
known equations, connecting deviators of stress and 
deformation s \- . f- and also their spherical components 

ct', e' , are fulfilled .s- = 2GT' , a' = K(3e'-BI ) , 
with the correction due to additional volume deformation, 
caused by neutron irradiation BI. Here G denotes shear 
modulus, K - volumetric deformation modulus. 

For those areas of solid body, where plastic deformations 
appeared the relation between deviators for simple loading 
can be represented as s'- = 2G')'-f'(f u , /, a ) ) . 

Here f'(e' u , I,a' k )~ the plasticity function, depending upon 
deformation intensity s' u , the value of neutron flow I and 

approximated parameters Cl' u . In the conditions of simple (by 

A. A. Iluyshin) loading [1] this function will be universal, i. e. 
it can be find from experiments with extension, torsion etc. 

So in the deformable body the relation between stresses and 
deformations under active loading from natural state and 
under influence of neutron flow in general case can be 
represented as 

s' li =2G3' t f(CI,a' l ), 

(3) 

& = K(3e' - B I) 

At that the plasticity function should be taken as 
I, Ct! k )= 1 in those areas, where s' u ^ s' s , s' s - 

deformation, corresponding to the plastic limit at the start 
time. 

If the force loading is rather quick (instantaneous) the 
irradiation reinforcement will not occur and originated areas 
of plastic deformations will be the same as in conditions 


without influence of neutron flow. Though if the active 
loading will be slow enough, external layers of the body will 
turn out to be reinforced and within these layers areas of 
plastic deformations will turn out to be smaller or will be 
missed at all, compared with the not irradiated body. There 
can take a place an effect when first plastic deformations will 
appear not on external reinforced surface, but under it, where 
the deformation intensity is great and plastic limit didn’t 
increased. 

So the influence of irradiation upon elastoplastic body is 
contrariwise to thermal, which decreases the plastic limit and 
results in increasing of areas of plastic deformations under 
equal loading. 

Let’s add to the relations in (3) differential equations and 
boundary conditions and also Cauchy proportions on the 
assumption of infinitesimal deformations 

o' ij , j +pF;=0\al j l j =R;onS (T , 

(4) 

u'i = u' 0i on S u ; 2s'- = u\ j +u' ji 

Comma in inferior index depicts the differentiation along 
next coordinate. Consider, that time variation of external 
loading and boundary moving occur in such a way, that 
appropriate loading trajectories are not related to the class of 
essentially complex loadings, and irradiation reinforcement 
takes place after force deformation of solid body. Hereinafter 
we shall suppose, that boundary problem (3) and (4) is solved. 

C. Problem definition for repeated alternating-sign 

loading 

Suppose that since t = 1 1 the influence of neutron flow 
disappears ((p = 0 ), and external forces change so that in all 
points of plastically deformable areas of the body V' takes 
place the unloading and following alternating-sign loading by 
volumetric F" and surface forces (at S a ) with 

ff ri 

boundary moving U gi (at o u ). The level of irradiation of the 
body remains constant and is equal to the value before 
unloading I l = (pt { . The plastic limit in points of the body 

depends upon coordinate z and becomes equal (z)), 

i. e. it depends upon the value of deformation and irradiation 
reinforcement. The scheme 

of the process discussed is shown at the fig. 2. 



Fig. 2 The scheme of the process of the deformation and 
radiation hardening 
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Let’s depict appropriate stresses, deformations and moving 
as O’": , £■' , u" . Formula (4) is valid for these values 

y y • 


c l » j + pF"~ 0; S^, 

u:=ul on S.M = 


(5) 


Relation between stresses and deformations will be n 
curve of the second cycle. 

s y =2G3 ii a" = 3^. (6) 

Here /"(e" , ej/, a" k )- plasticity function under repeated 
alternating-sign loading, depending upon the deformation 
intensity , preceding values of deformation intensity e\ 
and the level of irradiation of the body 1 1 , approximated 
parameters Cl'l , describing the deformation. 

At that plasticity function f is supposed to be equal 1 in 
those areas, where new plastic deformations didn’t appear, i. 
e. S u < £ v modulo, S - deformation, corresponding to the 

plastic limit <j" under repeated loading. 

Equations (5) and (6) define a boundary problem for values 
with two dashes. Complexity of the problem consists in the 
dependence of desired decision upon the unloading points 

( e\ , <j [ ), as the boundary problem must be defined and 
solved at every point the solid body. Let’s discuss one of 
methods to avoid such difficulties. 

For the values before the beginning of the unloading we 

shall retain designations CT-. , £.. , U j . Following Moskvitin 
[1] we shall define the following differences for moment 
t>t l : 


* 






(7) 


Let’s depict physical state equations for values with 
asterisks. In zones V' e and V" of unloading and elastic 
deformation are valid relations: 


4 = 2G3 i /-’ £ ! - e y(7i), 

in the area V" , where during variable loading the plastic 
deformation changes, must be valid relations: 

4 =2G 4 (8) 

Here / (e„ ,8 l ,I l ,a k ) , in general said, appears some new 

universal function, depicting nonlinearly of deformation 

* * 

diagram in axes <7 ~ £ (see the fig. 2). On the linear 
section / = 0 . 


o' = K(3e' - BI { ) , o" = K( 3s - Bl ) , that’s why for values 
with asterisks. 


CT =3Ke . (9) 

Equilibrium equations, boundary conditions and Cauchy 

* * * 

proportions for values CT , £ jf , M will be 

oI,j+pf*= o, f;=f;-f i " 

o } ;i i =R*',R; = R;-R:,onS a ; ( 10 ) 

4 = Ki = u' 0i - u " n , on S u ; 2 £* = u* + u] .. 


Relations (8)— (10) make a new boundary problem for 
values with asterisks. If we now suppose that the function / 
in all points of deformation curve can be approximated by 

function f , that is to depict it by the same analytical 

* 

equation but with another parameters a k , we shall exclude 

r * 

the dependence of / upon £ { : 

f*=f'( £ u’ 1 o<h)- 

Comparing after this equation, the formulas (3) and (4) for 
the body with loading from natural state and equations for 
values with asterisks (8) - (10) we can mark, that they 
coincide with closeness within designations. That’s why the 
solving of the problem for the values with asterisks can be get 
from known solving, of the problem, appropriated to loading 
from natural state by some replacements. For example, if the 
moving u\ = u\ (x, s', , e' s , I, a' k ) is known, then appropriate 

moving u { =u i (x,£ u ,s s ,I l ,a k ) and under repeated 

alternating-sign loading is calculated from equation (7) 
ff f * 

Uj — Uj - Uj . Stresses and deformations are calculated by 
formulas of the same type. 


III. RESULTS 

Obtained result may be extended in case of any given 
n-cyclic loading (theorem about cyclic loadings of 
elastoplastic bodies in neutron flux). Suppose that under 

n-loading by external forces F" , R" with boundary moving 
4 , stresses o"j , deformations and moving u ” appear. 

At the same time equilibrium equations, boundary conditions 
and Cauchy equations must be valid: 


ol J+P F,"=0-,a;i i = R:onS cr , 


V J 


u" = Mo; on S u ; 2e~ =u" j +u n jJ . 


(11) 


defined the following differences: 


/ -\\n /n 

Ojj =(-l) (Ojj 


n - 1 


o”), s/ =(-!)" (8, 


n — 1 




In all points of the body the volumetric deformation 
remains elastic. Consequently before the beginning of 
unloading and for current state equalities are valid 


*n s i\n / n — 1 n \ 

Uj =(-l) (Uj -Uj). 

Then the equations (11) turns out to be valid for the values 
with asterisks also: 
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( 12 ) 


° 7,+PF? =0, 

F* n = -f; 1 ); 

Vy" l j = R*" °n S a ,u* n = u* 0 'l on S u ; 

R; n =(-ir\Rr l -R';x 
C = (-!)"«■ 1 -<■); 

2 <" = <■+«;■• 

Let’s accept that under any n-loading the relation between 
sphere components of stresses and deformation tensors 
remains elastic. Repeating the previous supposition about 
possibility of curves s'y ~ 3^ and .y,” ~ a-" by nonlinearly 
functions of identical analytic type 


-2G 3 yf\z u n , !„<$), 


(13) 


we shall conclude that the solving for the problem for values 
with asterisks (12) and (13) under any given n-loading maybe 
get from the problem, concerning the loading from the natural 
state. For example, if the moving is known, 
u\ =u' t (x, s' u ,e' t ,I, a' k ) , then the appropriate value with 

asterisk will be m*" = u ■ (x, s*" , £*" , /, , ) . 

After this the desired moving u " can be calculated from 
equation 

«f=«;-X(-i)V*- (i4) 

k=2 

Stresses and deformations are calculated by formulas of the 
same type (14). 

As an example the problem of cyclic irradiation-force 
winding of sandwich beam with one embed end (fig. 3). 
Unsymmetrical with respect to thickness sandwich beam with 
external bearing layers made of metal and incompressible 
with respect to thickness internal layer (filler) made of 
polymer is considered. For the description of the pack 
kinematics the hypothesis of broken normal line is accepted: 
in bearing layers Kirhgoff hypothesis is valid, in the filler 
normal line remains rectilinear without changing the length, 
but it turns on some additional angle v|/(x) . Bearing layers 
are accepted to be elastoplastic, the filler - elastic. Analytical 
solution of appropriate problem of theory of elasticity is 
depicted in [3]. Solving of the problem of small elastoplastic 
deformations under loading from natural state was found by 
method of elastic decisions. Aluminium alloy was used as a 
bearing layer, and Teflon was used as filler. 



z 


y 



Fig. 3 cyclic irradiation-force winding of sandwich beam with 
one embed end 


The solution of this problem for the theory of elasticity is 
known [4]-[ 6]: 


y/ M (x) = C[ n> sh(j3x) + C 3 (n) + chi fix) + 
sh( y 5x)J 1 g t ' ,) ch (fix)dx -ch(fi x 
| g w sh{/3 x)dx 

u (n) (x) = y 3 i// (n) + 
-a 4 L 2 \p-p ( ;- l) ) + 

+ ai ^{q- q T l) ) + ^C[ n) x 2 

+ Cj n) x + Cg n) , (15) 

w ( "\x) = — [a x \y/ (n) dx- 
«2 J 

-a 1 L } \p-p^) + 

+ a x L-}(q-q^) + \a l Cl n) x^ + 
o 

+ -c ( ; ) x 2 + c\ n) x+c ( 6 n) x. 

2 4 5 6 




g (n \x) =—r l hr > + r 2 (p-p ( r ) ) + 


a. 




(n-l)'i 


a , 


r,(J(4-9r l> )*+c, w ). 

Here N" = \b^G k \m t s[ k 'dz, 
5 k = l K 

3 *=‘ K 

c 

Q co = 2 b Q G 3 J ojj//dz. . 


(16) 


Where n is approximation number; , h"~ , q'£ are 
called «additional» external stresses, assumed to be zero at the 
first step and are further calculated from the results of 
preceding approximation. For this aim, the formulas similar to 
(16) are used in which all summands have a superscript «n 
- 1 »: 




b , 


co(n-l) 


h<n-» = _ {H ^n-l)_ Qa (n-l) ) ^ 


b o 


co(«-l) 


N 


co(n-l) 


-X* 


(*)(«- 1) _ 


k=l 


= 3 Z, oZj G ^® W ( £ 

*=1 h, 


C*X«— 1) 1) 


dz. 


M “(""D = 


(*K«- 1) _ 


k = 1 
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= 2b 0 1 G 3 ro (3) (e^ 3)( ' i_1) )\|/ ( " -1) rfe , 

—c 

2 

Coefficients otj.o^.ou.p ,Yi>Y 2 >Y 3 and linearly 

here 

integrates operators L^,L^,I7^ are determinates in [4]-[6]. 

Under the boundary conditions for edge-fixing of the plate, 
we obtain the following recurrent formulas for the integration 

constants C{"\ ...,Cg' n : 

c, w = -Vte -«r 


C*' 0 = - 

ch (/?/) 

j g (n) tih(/] x)dx\ x=l 

- p 

sh (pi) 

-\g (n) Mpx)dx\ x= ^ 


-Jg' n V/l(/?.K)*/x|, =/ ], 

p 

C^=^L-\ P - P r%„- 

a 2 


—^i(q — < zi B- 1 ) )L=/ -—c[ n) i. 


a. 


a. 






a , 


« 2 


c )= — ^ 1 (p- J pr%=o- 

a 2 

r/ ryJ 




a 2 

Oj j- (n-1KI O' 




1 /-'Ml 


a , 


ff. 


(17) 


c<"' = r±L; 1 (p- P r , )l v . 


-*4 r-1 


(n-l)'l 


a. 


/ a=0 


“^3 (?“ ?») *=<>• 

a 2 

Under the boundary conditions for edge-fixing of the plate, 
we obtain the following recurrent formulas for integration 
constants: 

C{ n) =- i (L 2 \q-q^-%^-L 2 \ q -qt 1 % =l ), 


r (n) _a 7_( d <h)| _r)(")| ) 1 T~ l (n >1 

*-4 ~~ V'l |a= 0 |*=// ^2 ^co l|.r=0 > 


(«) _ n ( "'l 


cr=D 


-D 


M I 


2 x =0 ^2 a -/ 


Cljl 

2a, 


k°i 


- D ^ I 

x=0 u \ \x=, 


•()- 




c«— i) > 


2a, 


C<'° = 


(n)| 


_ n(")| 

A=0 ^1 A=/ 


A 


(n) _ a i 


aJ 




la 


2a, 


-^I^(p-p^), 


D 


la 

(n) a l 

/a- 


V - ’ 


ah- 




aj 

6a 2 


L 2 (q~q, 


(n-l)N 


These constants remain form 


^Li(p-pr)- 

la 2 

r (n) r (n) r (n) r {n) 


as in 


(17). 

Appropriate mechanical properties of mediums are 
depicted in [2]. On figures (4) and (5) shear y/ and flexure 
W of sandwich beam, calculated with the help of different 
physical state equations, are shown. Curves with one dash 
corresponds to loading from natural state, with two dashes - 
repeated cyclic winding due to alternating-sign loading: V - 
solving of elastic problem: 2— instant elastoplastic without 
irradiation; 3’- elastoplastic winding of previously irradiated 
beam (/,=5*10 24 m~ 2 ). 




Fig. 5 flexure W of sandwich beam 
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Under combined influence of force loading and neutron 
flow during interval t x till the value /, deformation will 
correspond to the curve 2'. Consequent instant unloading and 
force alternating-sign loading with level of irradiation /, will 
cause the shear and flexure of the beam, shown by curves 2". 
If the beam under cyclic loading was irradiated beforehand, 
then deformation would correspond to the curve 5". 


IV. Conclusion 

If you are using Word, use either the Microsoft Equation 
Editor or the MathType add-on (http://www.mathtype.com) 
for equations in your paper (Insert I Object I Create New I 
Microsoft Equation or MathType Equation). “Float over text” 
should not be selected. 
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